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Abstract 

We investigate via Monte Carlo simulations the kinetically constrained 
Kob- Andersen lattice glass model showing that, contrary to current expec- 
tations, the relaxation process and the dynamical heterogeneities seems 
to be characterized by different time scales. Indeed, we found that the 
relaxation time is related to a reverse percolation transition, whereas the 
time of maximum heterogeneity is related to the spatial correlation be- 
tween particles. This investigation leads to a geometrical interpretation 
of the relaxation processes and of the different observed time scales. 

1 Introduction 

The main challenge in the field of supercooled liquids is the understanding of 
the rapid increase of relaxation time and viscosity as temperature decreases pQ . 
Experimental evidences and theoretical models predict the simultaneous emer- 
gence of Dynamical Heterogeneities (DHs) [2 , 3 , which play the role of critical 
fluctuation in ordinary critical phenomena. Indeed, these studies depict DHs as 
dynamically correlated clusters: by lowering the temperature, the size of such 
clusters grows and their cooperative rearrangement becomes more complicated. 
This results in a rapid increase of typical cluster life- time, which is in turn re- 
lated to the structural relaxation time, at a macroscopic level. Being defined 
as the volume integral of a four-point correlation functon, #4(7*, t), the dynam- 
ical susceptibility x^)i estimates the volume of these clusters. Thus, %4 is 
expected to grow as some power of the dynamical correlation length, £4, derived 
by g4(r,t) ; in particular x<±(t) ^ f° r compact clusters in d dimensions, 

whereas a smaller exponent is expected if clusters have a fractal structure. The 
time £*, where X4 W reaches its maximum, %|, is also the time where correlations 
start to decrease, and thus it is interpreted as an estimate of the typical cluster 
life-time. Then, if DHs and relaxation are strictly related, t* ex r is expected. 
However, there exist systems where xX 1S found to decrease on approaching the 
transition [5j El H] , as well as early studies suggesting that the time of maximal 
correlation between particles displacements does not scale with the relaxation 
time r [7 j. Accordingly, the relation between the dynamical susceptibility and 
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the relaxation process remains elusive. Moreover, in order to unveil the precise 
relation between %4 and £4 one needs to know the explicit form of #4(7*, t), which 
is not an easy experimental task. 

Here we address this problem via a numerical study of the Kob-Anderson 
kinetically constrained lattice gas model [8 (KA-model), where it is possible to 
obtain very accurate data for the dynamical correlation length. After shortly 
reviewing the numerical model (Sec. we illustrate the relaxation process 
(Sec. |3| and the DHs (Sec. [3|. In Sec. [5] we show that the emergence of 
DHs is well described by the diffusing defect picture, which correctly predicts 
the exponents characterizing DHs. In Sec. [6] we explain the relaxation process 
in terms of "pacman" percolation, a model originally introduced to explain 
degradation of a gel due to the action of enzymes [9]. This picture furnishes a 
geometrical interpretation of the relaxation process and of the different observed 
timescales. This paper is based on the elaboration of results previously published 

2 Model and simulation details 

The KA-model [8 is a kinetically constrained model [13], which consists in a 
cubic lattice of volume V = L 3 containing N particles. Periodic boundary con- 
ditions are imposed along the x, y and z directions. The global density p = N/V 
is the only control parameter and it plays the role of an inverse temperature. 
No interactions between particles are present apart from an hard-core repulsion, 
which prevents more than one particle to occupy the same lattice site: all mi- 
croscopic configurations where the particles occupy N among V available sites 
are allowed, isoenergetic and equiprobable. Thus, the model is characterized 
by a trivial energy landscape, where no thermodynamic transition is available 
and time translational invariance holds. The system evolves following dyamical 
rules based on kinetic constraints: a continous time stochastic process allows a 
particle to move in a near empty site if has less than m = 4 neighbors, and if it 
will also have less than m = 4 neighbors after the move. Previous studies have 
shown that this model reproduces many aspects of glass forming systems: the 
dynamics slows down on increasing the density, suggesting the existence of a 
transition of structural arrest at p^ a — 0.881 [8] [14j [15j [16] . Nevertheless, it has 
been demonstrated that in the thermodynamic limit the transition of dynamical 
arrest only occurs at p = 1 [17] . 

Our simulations of the KA-model span a set of density values ranging over 
up p = 0.87. We have investigated a wide range of system sizes from L = 8 up 
to L = 50. Data reported in the following concern a system size L = 30, where 
we have performed the largest statistics: for each density values, the results are 
averaged on at least 10 2 over up 10 4 runs. The runs at higher density values 
last up 10 8 Monte Carlo sweeps. 

3 Structural Relaxation 

Here we investigate the relaxation process focusing on the time evolution of 
the density of persistent particles p(t) = y Yli=i n i(t)i where rii(t) is known as 
'persistence': rii(t) = 1(0) if site i is (is not) persistently occupied by a particle 
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Figure 1: Persistent particles in a numerical simulations of the Kob-Andersen 
model at p = 0.85, at times t x = 3.5 10 5 , t 2 = 7.5 10 5 , t 3 = 1.6 10 6 , and 
U = 2.1 10 6 . From Ref. [IT]. 




Figure 2: Normalized density of persistent particles (p)/p for different values 
of the density, as indicated. For p < 0.85, (p)/p is well described by the von 
Schweidler law, (p) / p = fo — {t/r) b , with / = 1 and b ~ 0.3. At short times, 
the dynamical susceptibility grows as t p , with p ~ 0.61. Adapted from Ref. 



in the time interval [0,t] [18]. p(t) is related to the high wave vector limit of 
the intermediate self scattering function [19 and describes the relaxation of the 
system in a very direct way: initially p(0) = p, but as time proceeds, particles 
eventually move from their original positions and p(t) decreases. This process 
is represented in Fig. [l] which shows the persistent particles in the system at 
different times of the same run. 

Accordingly, the normalized average value of such a correlation function, 
, represents the dynamic order parameter for relaxation. Fig. 2 shows 
that, in a large time window and for p < 0.85, p(t) is well described by the von 



3 



Schweidler law: 



^ = fo - it/r)\ (1) 

with b ~ 0.3 and /o — 1> whereas for larger time a stretched exponential fit works 
better. Being the relaxation time r defined as Mill = \je, Fig. j^D shows r 
to diverge approaching the transition of structural arrest, r{p) oc {p^ a — p) _Ar , 
with A r ~ 4.7, consistent with the result of Ref. [14]. 



4 Dynamical Heterogeneities 

Fig. [T] show that as less and less persistent particles survives, spatial correla- 
tions between them emerge. These correlations are quantified by the dynamical 
susceptibility X4W> related to the fluctuations of p, 

X4 (t) = Y - (ip(t) 2 ) - (P(t)) 2 ) , (2) 

and to the volume integral of the spatial correlation function between persistent 
particles at time t, 

I v 

X4(*) = ^X>M) (3) 

where 



54(r,t) = (ni(t) nj (t)) - {ni(t)){ nj (t)), r = \i - j\. (4) 

The spatial decay of #4(7*, t) defines the dynamical correlation length ^(t). In 
the following we describe the behaviour of #4(7*, £), ^(t) and X4^)- 



4.1 Four-point correlation function 

In Fig. [3|i we show #4(7*, t) at different times as a function of r. It clarifies that 
the spatial extension of the correlations grows until an intermediate time and 
then decreases at larger time. Its value at r = 0, #4(0, t) = (p(t))(l — (p(t)}), is 
determined by Eq. [4j 

The behaviour of #4(7*, t) is consistent with a scaling form such as: 

e r/UW 

^(r,t)=A(t) ;3= ^, (5) 

where is the dynamical correlation length (whose behaviour will be dis- 

cussed in the next session) and A(t) is the amplitude. 

Fig. ^jp shows #4(7*, t) normalized by its value at r — as a function of r, 
for different densities and at intermediate time (roughly of the order of t*). As 
expected, #4 becomes increasingly long ranged when the transition of structural 
arrest is approached. 
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Figure 3: Panel a: Correlation function between persistent particles #4(7*, t) 
as a function of the distance r at p = 0.85 and for different times t± = 10 2 , 
t 2 = 5.3 • 10 3 , tS = t* = 5.6 • 10 5 , U = 2.8 • 10 6 , t 5 = 4.6 • 10 6 . Panel b: 
#4(r, £ — as a function of the distance r at different values of the density. 
From Ref. [10]. 

4.2 Dynamical correlation length 

From the data relative to #4(7*, t), we have extracted the correlation length £4^) 
via an exponential fit of the initial decay. The behaviour of is illustrated 
in Fig. [4^i for different values of the density, and is well described by 

&(t) oct a exp(-at/tl) . (6) 

Accordingly, at short times £(£) grows as t a with a ~ 0.156 , and then it 
decreases after reaching its maximum value £| at time t| . We find that the time 

diverges as t| oc (p^ a — p) •> with At* = 3.8 ±0.1. 

From this result and from Eq. [6] it follows that the maximum dynamic cor- 
relation length diverges as £| oc t^ a oc (pk a — p)~ v with v = a\ t * ~ 0.54, in 
very good agreement with data (Fig. |5j. Note here a very important result: 
contrary to what is expected, we find that r and are not proportional but 
seems to diverge with different exponents, A r > \ t * , approaching the transition 
of structural arrest. The presence of a growing correlation length suggests that 
the system is approaching a critical point as the density increases. This sce- 
nario is conveniently described interpreting p = — log(t) as a chemical potential 
for the persistent particles, considering that the density of persistent particles 
monotonically decreases as time advances. The line where the correlation length 
reaches its maximum value in the p-p plane can therefore be interpreted as a 
Widom line, which in a second order transition ends at the critical point. The 
results of Fig. [6] suggest the presence of a critical point located at p = pka and 
p = — 00, where the correlation length diverges. Actually, the Widom line will 
possibly bend, ending at p = 1 where the transition is known to occur in the 
thermodynamic limit. Such an approach may open the way to a renormalization 
group treatment of the glass transition. 
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Figure 4: Panel a: dynamical correlation length for different values of the 
density. Panel b: divergence of the relaxation time r, of the time where the 
correlation length acquires its maximum value t|, and of the time where the 
dynamical susceptibility acquires its maximum value, t* . At low density, t* ex r, 
while at high density t* oc t|. Errors on t| and t* are of the order of 5%. From 
Ref. HI]. 
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Figure 5: Dynamical correlation length at t = t*, and prediction of the diffusing 
defect picture, ex t* a (xr q ,q = a\ r /\ t * . The full line is a (pk&-p)~ u , v — 0.54 
(we fix pka — 0.881 as estimated from the divergence of the relaxation time). 
From Ref. [IT]. 

4.3 Dynamical Susceptibility 

The emergence of an increasingly heterogeneous dynamics is clearly signaled by 
the dynamical susceptibility X4W? shown in Fig. [7] for different values of the 
density. Qualitatively the behaviour of x±{P) appears similar to the one observed 
for £4(t), initially growing as X4(t) cx t p , with p c± 0.6 and then decreasing after 
reaching its maximum value xt a ^ a time t* . The decoupling between t| and r 
strongly influences xj(p) and the time t* (p), leading to a complex behaviour. 
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Figure 6: Main panel: Widom line in the density, chemical potential plane. 
Circles indicates the time where the correlation length reaches its maximum 
value, at each value of the density. The continuous line corresponds to (1— Pka) oc 
(t/t*) a / l/ , and suggests that the system approaches a critical point at p — — oo, 
and p = pka- Inset: scaling of the dynamical susceptibility for 7 values of the 
density, in the range 0.78-0.87. From Ref. [TT] . 

In facts, inspired by Eq. [3] and Eq. [5] we find that x^(f) is well approximated 
by 

X,{t) oc Aim) 2 - 71 = [(P(*)>(1 - <p(t)»] ^i) 2 -", (7) 

where we have estimated A(t) = #4(0, t) and 77 c± consistent with Ref. [21]. Ac- 
cordingly, the behaviour of the susceptibility is essentially given by the product 
of two competitive factors, the amplitude and the correlation length. Indeed, 
at low densities t| ^> r, the amplitude dominates and Eq. [7| predicts t* oc r. 
By contrast, at high density t| <C r, and the maximum 01 the susceptibility 
X\ occurs at t* oc t| . Such behavior of t* is apparent in Fig. [4j at intermedi- 
ate density t* shows a clear crossover separating the two asymptotic regimes, 
where t* scales as r in the limit of low-density and as tg in the limit of high 
density. Asimptotically, when t* oc r, the maximum of the susceptibility scales 
as x2 oc r 2a oc (pfc a — p) -7 , with 7 = 2aA r , in agreement with our results. 
Conversely, when t* oc t|, we have X4 oc (tj) 2a oc (^|) 2 oc (p fca - p)~ q , with 
<7 2aA^*. 

5 Diffusing defects 

The results described so far are rationalized in the diffusing defects paradigm [2, 
[2TJ[22 , 23 , where the relaxation is ascribed to the presence of possibly extended 
diffusing defects, with density pd- The diffusing defects picture provides that at 
short times the correlation length grow as oc t x l x . The number of distinct 
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Figure 7: Dynamical Susceptibility, X4(t), for different values of the density, as 
indicated. At short times, X4(t) grows as t p , with p ~ 0.61. Adapted from Ref. 



sites visited by a defect increases as n v (t) oc t df / x , where df and 1/x are respec- 
tively the fractal dimension and the diffusion exponent characterizing the defect 
walk. Thus, the total number of distinct sites visited by the defects at time t 
is proportional to pdn v (t). These assumptions are expected to hold at least at 
short times, i.e. before defects interact. Under this condition, each site visited 
by a defect corresponds to a particle which first moves from its original position. 
Thus, we expect that at short times the total number of distinct visited sites 
scales as the density of particles which have already relaxed at that time (the 
non-persitent particles), 1 — M^-. From the decay observed for the density of 

persistent particles, we can infer 1 — (p(t)) / p = (l/r) b t b oc Pdn v (t) oc pdt df ^ x . 
Therefore this picture reproduces the von Schweidler law, and relates the density 
of defects pd with the relaxation time [23] . pd oc r~ b . By comparing the correla- 
tion length scaling provided by diffusing defect picture with the observed short 
time behaviour, ^(t) oc t a (Fig. [5|, we can conclude 1/x = a and df — b/a c± 2. 
Then we find that defects have a sub-diffusive nature, although they conserve 
the same fractal dimension of usual random walkers. We suggests this may 
be ascribed to defects that behave as random walkers, although characterized 
by a fat-tail waiting time distribution and, possibly, by weak spatial correla- 
tions which slows down the diffusion. In fact, it has been proved that these 
factors do not affect the fractal dimensions of the walkers [25j [26] . Dealing 
with the susceptibility, the diffusing defect picture predicts that at short times, 
X4(t) pdn v (i) 2 oc t 26 , which compared with our result allows to correctly 
estimate p = 2b. 



6 Reverse percolation 

Looking for a geometrical interpretation of the relaxation process, we were in- 
spired by the analogy with chemical gels. The mechanical rigidity of chemical 
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gels, in fact, arises from a percolating network of polymers which acts as a 
backbone in a liquid media. Since the bond of such network are of a chem- 
ical nature, relaxation cannot spontaneously occur and the system is solid at 
any time-scale. However relaxation may be induced by external factors, such 
as diffusing enzymes able to cut the bonds they meet [9j|27]: in this case, the 
system stays solid until the network survives, but it relaxes and becomes liquid 
at larger time-scales, when the network disappears. In our case the glass former 
may be though as rigid on time-scales smaller than the relaxation time, r. We 
suppose that during this time a percolating cluster of persistent particles plays 
the role of the physical backbone in gels; here the bonds are of a dynamical 
kind, i.e. we consider that two particles i and j are bonded in the interval [0, t] 
if they are nearest neighbours and persistent in this time interval. Defects, in- 
stead, play the role of enzymes and progressively destroy the cluster. A reverse 
dynamical percolation transition is expected for time-scales of the order of the 
relaxation time. In fact, as the absence of the percolating cluster leads to the 
loss of rigidity, we expect that this transition is related to the relaxation process. 

Data shown in Fig. [5] confirm our hipotesys. For density, beyond the onset 
of glassy dinamics, a cluster of persistent particles always spans the system at 
short times. We indicate with P(t) its strength, i.e. the density of persistent 
particles belonging to such cluster. (P(t)) vanishes at a time t per which is found 
to scale with the relaxation times, r, as the density increases (Fig. [8j inset). 

The figure also reveals that the cluster strength overlaps with the total den- 
sity of persistent particles, (P(t)) — (p(t)}, up to large times. This means that 
in this interval the percolating cluster is the only cluster present. At larger time 
p(t) slowly decays, while P(t) vanishes. This is due to the onset of finite clusters 
with a broad size distribution, that give contributions to p(t), but not to P(t). 
This circumstance may explain the crossover observed for the dynamic correla- 
tion function (see Fig. |2j: at short time, the decay of (p(t)) / p is characterized 
by a single relaxation time, which is the life-time of the percolating cluster, and 
this leads to the von Schweidler law (Eq. [I]). By contrast, the broad spectra of 
finite cluster life-time determine the stretched exponential decay at large times. 

To better understand the geometrical properties of this process, we investi- 
gate the correlation length, £ per (t), which is defined by the percolative correla- 
tion function g per (r,t): 

g P er(r, t) = g pc (r, t) - (P(t)) 2 = (n^n^t)) - (P(t)) 2 (8) 

where the pair- connected correlation function g pc (r,t) is limited to the pairs 
of particles that belong to the same dynamical cluster. g pc (r,t) can be also 
expressed as g pc {r,t) = P/-(r, t) + P-j(r, t) where P/.(r, t) and P£j(r,t) are 
the probabilities that two sites i and j belong to the same finite cluster, or 
to the percolating cluster respectively. Accordingly, if the percolating cluster 
is absent g pe r (r,t) = Pij(r,t) and f per (t) measures the typical size of finite 
clusters. Conversely, if finite clusters are negligible, then 

g per (r,t)=P°°(r,t)-(P(t)) 2 (9) 

measures the extension of the density fluctuations within the percolating cluster. 
In our case, at short time, (P(t)) — (p(t)} and P-j(r,t) ~ (ni(t)rij(t))\i_j\ =r , 
because almost all persistent particles belong to the percolating cluster, making 
the connectedness condition negligible [28]. Inserting these equalities in Eq. [9] 
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Figure 8: Percolation transition at p = 0.87. Left axis: dynamical correlation 
length £ (empty diamonds) and percolation correlation length £ per (full dia- 
monds). Right axis: density of persistent particles (p) (full line) and strength 
of the percolating cluster (P) (squares). The vertical dashed lines mark tg and 
tper, which is proportional to r (inset). From Ref. [TT] . 

and comparing it with the definiton of the four point correlation function, #4(7*, t) 
(see Eq. [4|, we find that g per (r,t) c± #4(7*, t), and consequently ^ per (t) — £4^)- 
Indeed, Fig. [5] confirms that the dynamical correlation length coincides with 
the percolative length, as long as finite clusters are negligible. The percolative 
length is affected by the two timescales characterizing the glassy dynamics, the 
time t = tg where the dynamical length reaches its maximum value, and the 
percolating time related to the relaxation of the system. At high densities, this 
makes £ per non monotonic, as in Fig. [5J 

7 Conclusion 

We have shown that in the KA model relaxation process and DHs are charac- 
terized by two different timescales r and tj, which implies that they are less 
tangled than expected. We explain this feature in the diffusing defect picture, 
where we relate the relaxation process to a reverse percolation transition, ob- 
taining a geometrical interpretation of the different time-scales. Accordingly, to 
their definitions, r occurs when a given large fraction of all sites has relaxed, 
while tg occurs when the correlations between the persistent particles decreases. 
Indeed, nothing forbids the correlations to decrease before a large fraction of all 
sites has relaxed, t| < r . However, in models of defects behaving as random 
walkers, one finds that r ex t| [2 . Therefore, one may speculate that, in our 
case, the decoupling between the two time-scales is due to a complex nature 
of defects. For instance, they may be non-conservative and characterized by 
birth and death rate with a constant average number, or they may not diffuse 
as random walkers. 
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